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ABSTRACT 

We develop an analytical method based on the lognormal approximation 
to compute the column density distribution of the Lyman-alpha forest in the 
low column density limit. We compute the column density distributions for six 
different cosmological models and found that the standard, C(95i?-normalized 
CDM model cannot fit the observations of the Lyman-alpha forest at z — 3. 
The amplitude of the fluctuations in that model has to be lowered by a factor 
of almost 3 to match observations. However, the currently viable cosmological 
models like the slightly tilted C05i?-normalized CDM+A model, the CHDM 
model with Vt v = 0.2, and the low-amplitude Standard CDM model are all in 
agreement with observations, to within the accuracy of our approximation, for 
the value of the cosmological baryon density at or higher than the old Standard 
Big Bang nucleosynthesis value of Q b h 2 = 0.0125 for the currently favored value 
of the ionizing radiation intensity. With the low value for the baryon density 
inferred by Hogan & Rugers (1996), the models can only marginally match 
observations. 

Subject headings: intergalactic medium - quasars: absorption lines - cosmology: 
large-scale structure of Universe 
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1. Introduction 

The Lyman-alpha forest, numerous weak hydrogen absorption lines in the spectra of 
distant quasars, have been the focus of extensive study for more that two decades. Recent 
cosmological hydrodynamic simulations greatly advance our understanding of it and allow 
unambiguous confrontation between observations with theoretical models (Cen et al. 1994; 
Miralda-Escude et al. 1995; Hernquist et al. 1996; Zhang et al. 1996). 

Numerical simulations show that the Lyman-alpha forest forms as a result of absorption 
of the quasar light by neutral hydrogen in inhomogeneities on small scales very much 
like the Lyman limit systems, arising from absorption by galaxies on larger scales. One 
of the predictions of such models is that the low column density Lyman-alpha forest 
(iVni ^ 10 14 cm -2 ) resides in regions of low overdensity and even underdense regions of 
the universe. While higher column density Lyman-alpha systems^ form in significantly 
overdense regions (5 = Sp/p ^ 3 — 10), which requires the use of numerical simulations 
to compute their properties accurately, the low column density systems may be studied 
analytically since there exist reliable approximations to study structure formation in the 
regime of low overdensity. The Lyman-alpha forest arising from these slightly overdense 
and underdense regions is the subject of this Letter. 

2. Column Density from the Stationary Phase Integration Around the Peak 

We first start with deriving an approximate expression for the column density of an 
absorption line arising from a density peak in velocity space. Let To be the Lyman-alpha 
optical depth for a homogeneous medium at the average cosmological density at redshift z 

3 We intentionally avoid using the term "clouds" since the forest mostly consists of 
fluctuations in the intergalactic medium rather than discrete objects. 



-4- 



along a line-of-sight to a distant quasar (see Jenkins & Ostriker 1991). We will consider 
the range of z ~ 3 when hydrogen was in ionization equilibrium with the radiation field 
which had photoionization rate V = 4.3 x 10~ 12 J_2i s _1 . J_2i, which characterizes the 
radiation intensity, is assumed to be 0.5. Since the number density of neutral hydrogen is 
proportional to the square of the density times the recombination coefficient, we can now 
compute the column density iVm of a Lyman-alpha absorption line: 



where £ = ln(l + 5), 5 being the overdensity, and the integral is taken over the spatial 
comoving coordinate x in the vicinity of the density peak. In this expression we implicitly 
assume that each Lyman-alpha line is dominated by a single density peak and we ignore 
peculiar velocity effects. The exact limits of integration are assumed to be unimportant 
as long as the dominant peak is included. The parameter a measures the deviation from 
the isothermality; if the recombination coefficient changes with the temperature as T^ 13 , 
and the equation of state is T ~ p 7 , then f2ni ~ p 2 T _/3 ~ p 2 ^ 13 " 7 and a = (3 r j/2. For 
typical values of (3 = 0.7 and 7 = 0.5 (as shown in Hui & Gnedin 1996), the correction for 
non- isothermality is only (3^/2 ^0.2. In the following, we will assume a = 0, but will allow 
for a range of temperatures, from 3, 000 K to 30, 000 K, at every value of density, to account 
for not only power-law equations of state, but also other relationships between the density 
and the temperature of the gas. However, we will keep a in our calculations for the sake of 
generality. 

The integral in (JJ) can be computed using the Stationary Phase method, 



where £" is the second derivative of the logarithm of the density with respect to the 
comoving coordinate along the line-of-sight and the subscript c denotes evaluation at the 
density peak. 




(1) 
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In order to proceed further and compute the column density distribution, we need to 
know the distribution function of the density and its derivatives. In the linear regime, the 
distribution function is gaussian, and the calculation can be performed fully analytically. 
However, we want to apply the approximation (||]) to a wider range of overdensities, 
particularly for underdense regions. We therefore use the lognormal distribution for the 
density (Bi, Borner, & Chu 1992; Coles, Melott, & Shandarin 1993), which is mainly based 
upon the assumption that the velocity distribution remains gaussian for larger values of 
overdensities than the density distribution does. Since the continuity equation in the 
expanding coordinates can be written as: 

dt adx v 1 ' 

where d/dt denotes the Lagrangian derivative and v is the peculiar velocity, we can see that 
the density distribution stays close to lognormal (i.e. £ is normally distributed) as long as v 
is gaussian. 

Since a one-dimensional slice of a three-dimensional gaussian random field is also a 
gaussian random field, we can use the BBKS formalism (Bardeen et al. 1986) to derive the 
column density distribution of the low column density Lyman-alpha forest. Let us introduce 
the quantity g as follows: 

9 = In (^=fif ) + *o(l - «) = 2(1 - a)(£ - £) - \ In - a)£"] , (4) 

where the parameters o"o, R*, and 7 (which appears below) are defined as in the BBKS 
paper, and £ = — cTq/2 is the average value of £. 

Then the column density distribution (number of absorption lines per unit column 
density per unit redshift) of the Lyman-alpha forest is given by the following simple 
expression: 

d 2 Af Ly - a c dn pk 

dN m dz N Kl H Jn (l +z) 3 d 9 ' () ' 



- 6- 



where dn^/dg is the one-dimensional comoving number density of peaks per unit interval 
of g, which can be easily computed using the BBKS methodology, and is given by the 
following expression: 

^ = \ r x e -Q(g, x) dx (6) 

d9 2(27r) 3 / 2 i?,a 3 7 2 v / 9/5 - 7 2 (1 - a) 

where 



o 



and 



A{g, x) = £ - e = 2(1 1 _ , (g + ~ ln((l - a)x 



The integral in (||) can be easily computed numerically for a given values of cr and 7. 

In order to test our approximation, we compute the exact column densities using 
the Zel'dovich approximation and Miralda et al. (1996) line identification algorithm and 
for each identified line we apply the Stationary Phase approximation. We find that the 
Stationary Phase column density is within a factor of two from and somewhat below (for 
higher column densities) the exact column density. 

Since our approximation assumes the lognormal density distribution, we test it against 
the Zel'dovich approximation, as shown in Fig.[l|. We show with open symbols the exact 
column density distribution (multiplied by the column density to better demonstrate the 
differences) from our Zel'dovich calculation (see Hui & Gnedin 1996 for details) together 
with the Stationary Phase calculation using the density distribution from the Zel'dovich 
approximation (filled triangles). [] The number of absorption lines drops drastically for 



4 The exact computation lies below the Stationary Phase + Zel'dovich approximation 
at low column density because the Miralda et al. (1996) line identification algorithm 
significantly underestimates the number of low column density lines (M. L. Norman, private 
communication) . 
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Fig. 1. — The column density distribution for the Standard CDM model with <t 8 = 0.7 
from the full calculation using the Zel'dovich approximation and Miralda et al. (1996) 
line identification algorithm (open circles) and the Stationary Phase calculation using the 
Zel'dovich approximation for the density distribution (filled triangles) versus approximate 
analytical treatment using the Stationary Phase calculation and the lognormal model for 
the density field (solid lines). The lower line corresponds to the same initial conditions as 
the Zel'dovich approximation (the same smoothing scale), and the upper line corresponds 
to the smaller smoothing scale that gives the same rms density fluctuation as the Zel'dovich 
approximation. The number of lines drops drastically for the simulated distributions at low 
column densities due to finite resolution. 
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iVm < 10 12 cm~ 2 due to the finite resolution of our simulations. The two solid lines show 
the approximate analytical calculation for two choices of the smoothing scale: the lower 
curve has the same smoothing as the Zel'dovich approximation, and the upper curve has a 
smaller smoothing scale so as to have the same rms density fluctuation as the Zel'dovich 
approximation. In general, we again see that the analytical approximation gives a good 
guess for the real column density distribution, perhaps, slightly underestimating the column 
density distribution. 

3. Testing the Models 

We consider six different models whose parameters are compiled in Table |l|. Our 
treatment only requires three spectral parameters, &q, 7, and R*. However, for a power 
spectrum which behaves like k~ 3 at large k, integrals for 7 and diverge. The power 
spectrum, therefore, has to be cut off at some scale. The natural cut-off scale would be 
the Jeans scale but the power spectrum at scales smaller than it behaves like k~ 7 (Bi et 
al. 1992), and integrals needed to compute R* and 7 diverge. We, therefore, apply the 
exponential cut-off at the Jeans scale, 

P g Uk) = PMk)e-( k / k rf (9) 

where 

kj = 7.4(fi (l + z)W* K/Tj) 1/2 h Mpc" 1 , (10) 

but we treat the temperature Tj at which we compute the Jeans scale as an independent 
parameter not necessarily equal to the gas temperature T gas . We note here that this cut-off 
of the power spectrum is equivalent to smoothing the one-dimensional density distribution 
on the scale of Doppler broadening. Even though Doppler broadening does not affect 
the value of the column density, the identification of a line depends on the value of b 
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(b = ^jlksT jm v where m p is the proton mass) as two narrow peaks that are separated 
by a distance much less than b in velocity space would not be counted as two lines but 
rather as one line. Since this procedure is somewhat dependent on the line identification 
algorithm, we use the freedom in Tj as a way to account for differences between different 
line identification algorithms. We then consider the range of temperatures from 3 x 10 3 K 
to 3 x 10 4 K for both Tj and the gas temperature T gas in our calculations and we put a = 0. 
The range of temperatures at a given density should be enough to cover possible variations 
for any other reasonable equation of state. We assume J_ 2 i = 0.5 throughout. 

The first two models we consider are standard CDM models, one is with a$ = 0.4 at 
z = and the other is CCRE-normalized. We use the BBKS transfer function to compute 
spectral parameters <To, 7, and R* for this model. 

The next three models are the COBE-normalized CDM+A models whose transfer 
functions cannot be fit by the BBKS formula accurately. We therefore compute transfer 
functions for LCDM models using the linear gravity code similar to (but different from) 
COSMICS package (Bertschinger 1995) and we use exact transfer functions to compute the 
spectral parameters. The first LCDM model assumes a small tilt in the primordial power 
spectrum and 25% of gravity waves added in quadrature as in Kofman, Gnedin, & Bahcall 
(1993). 

Finally, we study the COBE- normalized CHDM model with currently favored value of 
Q v = 0.2 and no tilt. We use the Ma (1996) transfer function in our calculation. Table |2] 
lists spectral parameters for the six models at z = 3 and T = 10 4 K. 

FigfJ shows the computed column density distributions for the above six models where 
both Tj and T gas independently take values from 3 x 10 3 K to 3 x 10 4 K at z — 3. We use 
the old Standard Big Bang Nucleosynthesis value of VL^h 2 = 0.0125 for the baryon density 
in these calculations. We emphasize again that since the temperature of the intergalactic 
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Table 1: Cosmological Models 



Model 


^0 






h 


n 


0"8 


SCDM1 


1 








0.50 


1 


0.40 


SCDM2 


1 








0.50 


1 


1.15 


LCDM1 


0.35 


0.65 





0.70 


0.96 


0.67 


LCDM2 


0.35 


0.65 





0.70 


1 


0.85 


LCDM3 


0.40 


0.60 





0.65 


1 


0.84 


CHDM1 


1 





0.2 


0.50 


1 


0.76 



Table 2: Spectral Parameters at Tj = 10 4 K 

Model a 7 R*{h~ l Mpc) 

SCDM1 1.01 0.52 0.105 

SCDM2 2.90 0.52 0.105 

LCDM1 1.11 0.51 0.180 

LCDM2 1.50 0.52 0.178 

LCDM3 1.51 0.52 0.167 

CHDM1 1.18 0.49 0.106 
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Fig. 2. — Comparison between six cosmological models and observations (solid squares) at 
z = 3. The shaded region shows the predicted range for the models with both T gas and Tj 
taking values between 3 x 10 3 K and 3 x 10 4 K. 
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gas in not known precisely, and it may also depend on local conditions, it is appropriate to 
consider the hatched £1X6 cl cLS the region where a more accurate calculation would lie. The 
width of the hatched area also corresponds to the accuracy of the analytical approximation 
as can be judged from Fig.|l| We also show the observed abundances from Hu et al. (1995) 
with filled squares. The size of the squares roughly corresponds to the observational 
uncertainties. 

As could be expected, the COBE-normalized standard CDM does not fit the data, but 
the low-amplitude SCDM, slightly tilted LCDM and the CHDM models provide a good 
fit to the data for low column densities. At higher column densities, iVm ^ 10 14 cm -2 , 
our approximation breaks down since those column densities correspond to significantly 
overdense regions. The two no-tilt LCDM models produce slightly too few low column 
density Lyman-alpha systems because their o"o is too high and their underdense regions are 
too underdense; they may be only marginally consistent with observations. 

It is remarkable that the column density distribution is by far most sensitive to the 
value of o"o- It is the value of o"o(Tj = 10 4 K) « 1.1 that ir responsible to the match with 
observations for SCDM1, LCDM1, and CHDM1 models. We found that the dependence of 
the results on other parameters (Qo, h, J-21, T, 7, and i?*) is weaker than dependence on 
(To, and for any set of those six parameters (within reasonable limits) it is possible to adjust 
the value of <j to give a good fit to observations. 

Finally, we address the question of uncertainty in the baryon density. Recently there 
have been new determinations of the cosmological baryon density using the high-redshift 
deuterium measurements giving two contradicting values of f^/i 2 = 0.006 (Hogan & Rugers 
1996) and Vt b h 2 = 0.024 (Tytler & Buries 1996). We plot the results for the three best-fit 
models, LCDM1, CHDM1, and SCDM1, for those two values of the baryon density and 
J_2i = 0.5 in Fig.|3|. We note that the high value for the baryon density improves the 



-13- 



Q h h 2 =0.006, J_ 21 = 0.5 



n b h 2 -0.024, J_ 21 = 0.5 



10 9 hi 



s io 



-10 



e 

i 

-a 



io- 
io- 

io- 13 y 





i i i i ii iij 


i 1 1 1 in i i 1 1 1 nij 




i 












■ 








M ■ 








! : ■ 








: ■ ■ 








: ■ 


■ - 




LCDM1, 


n = 0.96 




tnl 


i 


i i' 


1 



ttt 1 — I I I I I I II 1 — I I I I I I II 1 — I I I I I I 1 1 r 



LCDM1, n = 0.96 



rrrn — i — i i — i — i mini — i — i i i — i — q nm — i — i i — i — mm 



10" 



N 



53 10 



10 — 



X5 



o 

i 



6 

o 

N 

i— i 

x 

a 

i 

w 
T3 



10-n L- 
lO" 12 
IO" 13 



CHDM1, Q v =0.2 



l-l- 1 1 1 I I I I I I tLLLl I I I I I I I 3) 



J I 

H 1 I Mill 



CHDM1, Q = 0.2 



io- 9 M 

10 -io 

IO" 11 
IO" 12 

io- 13 



V,\[ 1 


i 1 1 1 in i i 


1 1 1 lllj 


1 1 


1 1 • lllll: 










: ■ 








; ■ 








: ■ 










, ■ 

■ 




I SCDM1, 


cj 8 =0.4 




■ : 


tnl i 


1 1 1 


1 


' 1 1 



mri r 



my 



i i i in r 



SCDM1, (7 R =0.4 



J I I 



I L 



IO 11 10 12 10 13 

N HI (cm-*) 



10 14 IO 11 



10 12 10 13 10 14 

N HI (cmi 



Fig. 3. — Comparison between three best-fit cosmological models and observations for two 
currently favorable values of the cosmological baryon density. 
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agreement between the models and the data, while with the low value for the baryon 
density, the models are only marginally consistent with observations. 

4. Conclusions 

We have developed an analytical method to compute the column density distribution 
for the Lyman-alpha forest in the low column density limit, Nn I ^ 10 14 cm -2 . We computed 
the column density distributions for six different cosmological models and found that the 
standard, COBE-normalized CDM model cannot fit the observations of the Lyman-alpha 
forest at z — 3. The amplitude of density fluctuations in that model has to be lowered by 
a factor of almost 3 to match observations. The CO BE -YioimsXized CDM+A model with 
a slight tilt, the CO-RE-normalized CHDM model with Q u = 0.2, and the low-amplitude 
Standard CDM model are all in agreement with observations for a range of values for 
the cosmological baryon density and for the currently favored value of the ionizing 
background radiation intensity of J_2i = 0.5. The no-tilt LCDM models without gravity 
wave contribution produces slightly too few low column density Lyman-alpha systems. 
We also note that since at high redshift, z ^ 2, all cosmological models have similar rate 
of growth of perturbations on scales of interest (including the CHDM model, since the 
neutrino free-streaming scale is larger than the characteristic scale for the low column 
density Lyman-alpha forest), and since the three models, LCDM1, CHDM1, and SCDM1, 
have similar predictions at z — 3, there will be little differences in their predictions at 
z = 2 or z = 4. We conclude that, based on the lognormal approximation, the column 
density distribution alone provides only a weak constraint on currently viable cosmological 
models unless the physical state of the intergalactic medium is known more precisely 
from independent observations. A more accurate approximation would also be useful in 
decreasing the width of the hatched area in Fig.|2|, which takes into account our ignorance 
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of the conditions of the intergalactic medium as well as the accuracy of the lognormal 
approximation. 

The computed column density distribution is by far most sensitive to the amplitude of 
fluctuations at the Jeans scale, a . We conclude that models with a ~ 1 (at Tj = 10 4 K 
and z = 3) are capable of reproducing the observational data for a large range of other 
parameters. 

Finally, we note that with the low value for baryon density as inferred by Hogan & 
Rugers (1996), all of the models tested here are only marginally consistent with observations. 

The authors thank Edmund Bertschinger for numerous fruitful discussions and valuable 
suggestions. We are also grateful to Sergei Shandarin for very helpful comments. This work 
was supported in part by the NSF grant AST-9318185. 
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